HADRON AND QUARK FORM FACTORS IN THE 
RELATIVISTIC HARMONIC OSCILLATOR MODEL 



V. V. Burov 

Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 
Head Post Office, Box 79, Moscow, U.S.S.R. 

A. De Pace 
Istituto Nazionale di Fisica Nucleare - sez. di Torino, 
Via P. Giuria 1, 1-10125, Torino, Italy. 

S. M. Dorkin 

Far-East State University, Sukhanova 8, 
Vladivostok, U.S.S.R. 

P. Saracco 

Istituto Nazionale di Fisica Nucleare - sez. di Genova, 
Via Dodecanneso 33, 1-161 46 Genova, Italy. 

Abstract 

Nucleon, pion and quark form factors are studied within the relativistic harmonic 
oscillator model including the quark spin. It is shown that the nucleon charge, magnetic 
and axial form factors and the pion charge form factor can be explained with one oscillator 
parameter if one accounts for the scaling rule and the size of the constituent quarks. 
PACS: 12.35H, 13.40F 

In this paper, we will analyze the proton electromagnetic and weak form factors (FFs) and 
the pion charge FF in the framework of the Relativistic Harmonic Oscillator Model (RHOM). 
Studies in this model were carried out in refs. |1T[|-[|TT||. 

In ref. a method was proposed for building a covariant and gauge-invariant current 
within the U{12) ® 0(3, 1) model. There remain still open problems in that approach (see, for 
instance, lectures [jl3[)- In refs. @, ||, 0| a covariant and gauge-invariant current was found 
for the SU{6) ® 0(3, 1) model and it was shown that all the nucleon FFs can be described in 
this case. However, the agreement with the experimental data is here much worse than within 
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the nonrelativistic model that takes account of the Lorentz contraction of the nucleon wave 



function |]I4[, although this latter model fails to describe the electric FF of the neutron. 

The aim of the present note is to describe the nucleon FFs within RHOM using the SU (6) ® 
0(3, 1) scheme of derivation of the covariant and gauge-invariant currents under the assumption 



that the behavior of the FFs when — > oo is governed by the quark counting rules ||T^ and 



that there holds the experimentally observable scaling, i.e. G^^{q^) = G\j{q^) / = G\j{q^) / 
where ^n,p are the magnetic moments of neutron and proton. 

Enforcing the scaling law leads one naturally to introduce FFs for the constituent quarks. 
Actually, these could be regarded as interpolating functions between the static quark model 
and the asymptotic perturbative regime, since, as we shall see, they go to 1 both a = and 
oo. 

Let us consider a system consisting of quarks in the field of a relativistic harmonic 
oscillator. The corresponding wave function can be represented in the form 

)(a;i, . . . , x^) = A^^ixu X2, . . . x^)t/(^)(P )• (1) 

where A is the operator of antisymmetrization of quarks including the color degrees of freedom 
to be not written for simplicity, $Ar(xi,X2, . . -x^) is a covariant space - time wave function, 
and U^'^^ (p ) is a spin wave function to be described below. Let the wave function $ at obey 



the Klein-Gordon equation with a relativistic harmonic oscillator potential Ipjl - |11 




N N-1 

^ ^ (X, - XjY 

i>j j=l 



$iv(xi,X2, . ..Xn) = 0, (2) 



where Pi = —id/dxi is a 4-momentum, /C is the oscillator parameter, Xi is the 4-coordinate of 
the 2-th quark (we assume all quark masses equal because of isospin invariance). Passing to 
the center-of-mass coordinates X and the internal variables tq, . . . r^-i and diagonalizing, one 
can represent (0) in the form 

(p^ - Mp)$Arg(ro,ri, . . .rAr_i,p) = 0, = -2aAra+ai^ + const, = JCNVn, (3) 

where a^^ and af^ are, respectively, particle creation and annihilation operators. Under the 
Takabayashi condition 0, necessary for removing nonphysical oscillations along the coordinate 
of relative time, p^af^^Nq = 0, one gets the following solution 

$iv,(ro,ri,...r^_i,p) = (^)''"' exp(^i^^'^ ^V,,r,.), (4) 
^n{xi,X2, . . .Xn) = exp[ip^X^]^Nq{ro,ri, . . .rN-i,p), (5) 



2 



where p is the total momentum of the system and K^^'^ = g^^^ — l-p^-p^ j-p^ . 

The spin wave function can be constructed in two ways (see refs. [|], |10|) by transforming 
the nonrelativistic spin wave function in the rest frame. The first method is to transform 
the wave function of every quark separately as a Dirac spinor on the basis of the Bargmann- 
Wigner equation. The second is to transform the wave function of the system as a whole, the 
minimal Pauli transformation, since in this case the wave function has a minimum number of 
components. The latter seems favorable, as it allows for a non-identically-zero neutron charge 
FF, in contrast to the first method. So, we will follow the second method according to ref. [ p^ . 
Then, the spin wave function \J^^^ (p ) can be represented in the form 



f/W(p) = 5(p)f/(^)(0), f/W(0)=(^ 



(6) 



exp 



2|p 



exp \pxhH\ , 



Pi 




where x is the nonrelativistic spin function of the system, if=(p-CT)/2|p|,6 = cosh^^po/^ 
and a = Y^iLi ^\ being the Pauli matrices of the i-th quark. 

Based on refs.0, |^, 0, we write the electromagnetic action in the form 

Iem = J l[dx^T.Mxi, . . . ,XN)A^{Xk) = j dX j\^\x) A^{X) , (7) 



i=l 



where 



.Xn) 



d 



dx 



k 



d 



dx 



+ 



d 



kv 



dx 



kv 



p ■ 



In eq. (||) (^p^"*) is the initial (final) wave function of the quark system (p, Ck is 

the charge of the k-th quark, cr^^ are the spin matrices of the k-th quark (cr*j- = £y7crf , 
(7^4 = crfj = pidf ). Assuming the constituent quarks to be not point particles, we identify gsiq'^) 
and gKiiq^) with charge and magnetic quark FFs (g = p' — p is the 4-momentum transferred for 
an N quark system with an initial (final) 4-momentum {p'^))- Inserting the wave function (|^) 
into eqs. (0) and (H) and computing the integrals over the internal quark variables tq, . . . , ttv-i, 
one derives the matrix elements of the effective current for an N quark system between states 
of momentum p^ {p'^) and spin component s (s'): 



p'sVir\o)\ps) = ^^^E(f/f )(p')r.,,t/f )(p)). 



(9) 
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where 



(10) 



Here the overlapping integrals over space-time variables are the following 



lN{q^ 



1 



(1 + qV2M%) 
1 + Nq^/2Ml 
l + gV2M2, ' 



(Ar_i) exp 



A^- 1 



q' 



N , 



fill 



[12) 



Mn being the mass of the A^-quark system. Using (|^) one can write down eq. (|^) for a three- 
quark system, finally obtaining expressions for the nucleon FFs, namely 



F^{q'] 



= GW) 



9M{q')I^'\q' 



l + ^\h{q')9E{q')-^,9M{q') 



2M2 



4M2 



9M{q')I^'\q' 



4M2- 



I^'\q% 



(13) 
(14) 
(15) 



It is also possible to compute the nucleon axial FF FA{q^) = G'a(q'^)/G'a(0) (for details see 
ref. 0), 



FA{q') = GA{q')/GA{0) = 9A{q'W)I^'\q'), 



(16) 



where 9A{,q^) is the constituent quark axial FF. 

Note that within the SU{Q) model the magnetic moment of the proton equals yUp = 3, 
whereas for the neutron /i„ = —2, somewhat different from the experimental values /Xp = 2.793 
and /i„ = —1.913. Therefore, in the analysis of the experimental nucleon FFs we will normalize 
them to 1 at g2 = Q. In expressions (|l^ - |16D M = M3 = 0.938 GeV is the nucleon mass. 

Let us go back to the magnetic and charge nucleon FFs. According to the quark counting 
rules [0 the nucleon FFs decrease as when qq. From inspection of formulae (p!ll-|T6|), 
it is clear that all the FFs calculated in this model have the correct asymptotic behaviour 
already with point-like constituent quarks. For instance, when 00, Flji^q^) ~ 9M{q'^)/q^, 

therefore it is natural to make the minimal assumption that 9M{q^) = 1- However, since from 
the scaling condition it also follows that F^j{q'^) = F^lq"^), we can easily derive the following 
expression for the quark charge FF 9E{q'^)'- 

(1 + 3gV4M2) (1 + g72M2^ 



9E{q 



;i + gV4M2) (1 + 3g72M2 



:9A4{q 



(17) 
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We stress that eq. {\17\} is dictated by the scahng condition, which would be violated for point- 
like quarks, and not by the asymptotic dependence on q of the FFs, which comes out naturally 
from the model. 



The only free parameter in this model is the oscillator parameter: fitting (13) to the exper- 
imental data for the proton FFs, we find as = 0.42 (GeV/c)^, in agreement with the results of 
ref . . From formula (0) one determines the theoretical value of the slope of the electric FF 
at = to be —dF^{q^)/dq^\q2^Q = 0.022, which is consistent with the known experimental 
value 0.0202 ± 0.0003 fm^ [|15l. 

The standard parametrization of the nucleon FFs is in terms of a dipole function 

FD{q') = - ^ ^, 

1 + gVA| . 



where for the nucleon electromagnetic FFs A|; = 0.71 (GeV/c)^ and for the axial nucleon FF 
Aa = (1.032 ± 0.036) GeV/c |§. 

In Figs. 1 we show the comparison of our computations within RHOM with the experimental 
data for the nucleon electromagnetic FFs divided (except for Fg) by the dipole parametrization 
([T8|). Note that the proton FFs have been fitted (adjusting 03) to experiment, whereas the 
others have been found without additional parameters. In the case of the neutron, the standard 



Galster parametrization [|1^ is also shown: actually, by using an effective mass of ~ 1 GeV in 
the RHOM calculation one can closely reproduce the Galster curve. 

In Fig. 2 the dipole parametrization (solid lines) of the axial FF is compared to formula 
(17) with QAi'f) = 1 (dotted line) and QAi'f) = 9E{,<f') (dashed line). Note the improved 
agreement in the latter case: the effect of the constituent quark size, although not dramatic, is 
not negligible (~ 10% ^ 15%). 

It is straightforward to estimate the nucleon and quark sizes from equations (p!3D, (^4|) and 



^{P,n) = V< ip,n) > = V6(V2a3 + 1/M2), f, = ^< > = y'3/M2. Thus, f(p,„) = 0.74 
fm ('^(^p^) — (0.87 ± 0.07) fm), and fq = 0.36 fm. From (|15]), on the other hand, one finds 
< r2 >E^„= -3/M2 = -0.13 fm2 (< >|^„= (-0.119 ± 0.004) fm^). 

Concerning the axial radius, it is interesting to note that in the qa = 1 case V< >a = 



3/03 = 0.53 fm, whereas for ^^(g^) = 9E{q'^) V< r'^ >A = V3(l/a3 + = 0.64 fm (to 

be compared with the experimental value of (0.65 ± 0.07) fm). 

The procedure outlined above can also be applied to the quark-antiquark system. Here a 
source of ambiguity is due to the two possible representations for the booster in eq. (H): indeed, 
both a = and a = — a^ lead to a valid formulation 0. However, the latter case 

is known to give a simple pole-like asymptotic behaviour for the pion FF only for a special 
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combination of the single-quark currents in eq. @ (see ref. 0]). In the following we shall stick 
to the first formulation, where the pion charge FF reads 



FM') = 9E{q'W)I^'\q'). (19) 



gE^ff) being still given by eq. ([T7|) . 

Assuming the same oscillator parameter /C for the two- and three- quark systems, the value 
of ^2 can be derived from as since, using eq. (^), a2/«3 = (2/3) "^^/^^ Then, from 03 = 0.42 
(GeV/c)^, one gets 0:2 = 0.23 (GeV/c)^. Note that in the literature both 02 and 0:3 have usually 
been regarded as free parameters (see, e. g., |^). A known problem is connected to the mass 
of the quark-antiquark system: indeed, using the physical pion mass one largely underestimates 
F7r(g^). In the literature a wide range of "^[/(G) symmetric masses" has been employed (0.5 
GeV< M2 < 0.8 GeV) @, 0, |, § and in the following we use M2 = 0.77 GeV. However, it is 
found that this is not a critical parameter, since expression (T^) is mildly dependent on M2 for 



0.4 GeV < M2 < 1 GeV, giving in this range descriptions of the data of comparable quality 
(note that the same is not true for the case with qe = 1, which is much more sensitive to M2). 

In Fig. 3 we compare with the experiment our calculations of Ft^^q^) using either gE{<f) given 
by eq. (0) and gE^^f) = 1- The improvement due to the constituent quark FF is apparent. 
It is again interesting to note that the pion radius is independent of M2 and turns out to be 



= V< >n = y^3(l/2a2 + = 0.62 fm and = Y'3/2a2 = 0.50 fm for gE{q^) given 

by eq. ([T7|) and gE = 1, respectively (f^^^ = (0.66 ± 0.01) fm). 

We summarize here our conclusions. This model gives a simple description of the exper- 
imental data on nucleon and pion FFs provided only one arbitrary parameter is used. Note 
that the quality of the description can be improved upon by using guiq^) and gA{,q^) to fit the 
nucleon FFs. The model could then be used, with all the parameters fixed, in investigations, 
for instance, of the strong vrNN vertex, A - isobar FFs, nucleon structure functions and so on. 
Besides, the model can be easily extended to systems with a number of quarks larger than 
three, e. g. to the study of the deuteron FFs |11[]. 
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Figure captions 



Fig.l The solid line is the ratio of the normalized electromagnetic FFs of the nucleon to the 
dipole parametrization (p!8D; in the case of the neutron the squared electric FF is shown 
(solid line), compared to the Galster parametrization = — /i„(g^/4M^)FD (dashed 



line). The experimental data are taken from refs. ||17|| . 



Fig. 2 The normalized axial FF of the nucleon F^iq'^) = Ga{(1^)/G{()). The solid lines represent 
the dipole fits of the neutrino experiments and correspond to the upper and lower world 
average values for []TB| ; they are compared to the present calculation in the framework 
of RHOM with gA{q^) = 1 (dotted line) and gA^q^) = gsiq^) (dashed line). 

Fig. 3 The pion charge FF (^) with gsiq^) either given by eq. ([T7| ) (solid line) or equal to 1 
(dashed line). The experimental data are taken from refs. [Il9| . 
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Fig. 1 
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